Subsurface flow measurements using passive flux meters in variably-saturated cold-regions landscapes

To date, passive flux meters have predominantly been applied in temperate environments for tracking the movement of contaminants in groundwater. This study applies these instruments to reduce uncertainty in (typically instantaneous) flux measurements made in a low-gradient, wetland dominated, discontinuous permafrost environment. This method supports improved estimation of unsaturated and over-winter subsurface flows which are very difficult to quantify using hydraulic gradient-based approaches. Improved subsurface flow estimates can play a key role in understanding the water budget of this landscape.     

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.     

Geostatistical features of streambed vertical hydraulic conductivities in Frenchman Creek Watershed in Western Nebraska

This study evaluated the spatial variability of streambed vertical hydraulic conductivity (K_v) in different stream morphologies in the Frenchman Creek Watershed, Western Nebraska, using different variogram models. Streambed K_v values were determined in situ using permeameter tests at 10 sites in Frenchman, Stinking Water and Spring Creeks during the dry season at baseflow conditions. Measurements were taken both in straight and meandering stream channels during a 5 day period at similar flow conditions. Each test site comprised of at least three transects and each transect comprised of at least three K_v measurements. Linear, Gaussian, exponential and spherical variogram models were used with Kriging gridding method for the 10 sites. As a goodness-of-fit statistic for the variogram models, cross-validation results showed differences in the median absolute deviation and the standard deviation of the cross-validation residuals. Results show that using the geometric means of the 10 sites for gridding performs better than using either all the K_v values from the 93 permeameter tests or 10 K_v values from the middle transects and centre permeameters. Incorporating both the spatial variability and the uncertainty involved in the measurement at a reach segment can yield more accurate grid results that can be useful in calibrating K_v at watershed or sub-watershed scales in distributed hydrological models.     

The stable isotope signature of kilometre-scale fracturedominated metamorphic fluid pathways, Mary Kathleen, Australia

Abstract Large calcite veins and pods in the Proterozoic Corella Formation of the Mount Isa Inlier provide evidence for kilometre-scale fluid transport during amphibolite facies metamorphism. These 10- to 100-m-scale podiform veins and their surrounding alteration zones have similar oxygen and carbon isotopic ratios throughout the 200 × 10-km Mary Kathleen Fold Belt, despite the isotopic heterogeneity of the surrounding wallrocks. The fluids that formed the pods and veins were not in isotopic equilibrium with the immediately adjacent rocks. The pods have δ¹³C_calcite values of –2 to –7% and δ¹⁸O_calcite values of 10.5 to 12.5%. Away from the pods, metadolerite wallrocks have δ¹⁸O_whole-rock values of 3.5 to 7%. and unaltered banded calc-silicate and marble wallrocks have δ¹³C_calcite of –1.6 to –0.6%, and δ¹⁸O_calcite of 18 to 21%. In the alteration zones adjacent to the pods, the δ¹⁸O values of both metadolerite and calc-silicate rocks approach those of the pods. Large calcite pods hosted entirely in calc-silicates show little difference in isotopic composition from pods hosted entirely in metadolerite. Thus, 100- to 500-m-scale isotopic exchange with the surrounding metadolerites and calc-silicates does not explain the observation that the δ¹⁸O values of the pods are intermediate between these two rock types. Pods hosted in felsic metavolcanics and metasiltstones are also isotopically indistinguishable from those hosted in the dominant metadolerites and calc-silicates. These data suggest the veins are the product of infiltration of isotopically homogeneous fluids that were not derived from within the Corella Formation at the presently exposed crustal level, although some of the spread in the data may be due to a relatively small contribution from devolatilization reactions in the calc-silicates, or thermal fluctuations attending deformation and metamorphism. The overall L-shaped trend of the data on plots of δ¹³C vs. δ¹⁸O is most consistent with mixing of large volumes of externally derived fluids with small volumes of locally derived fluid produced by devolatilization of calc-silicate rocks. Localization of the vein systems in dilatant sites around metadolerite/calc-silicate boundaries indicates a strong structural control on fluid flow, and the stable isotope data suggest fluid migration must have occurred at scales greater than at least 1 km. The ultimate source for the external fluid is uncertain, but is probably fluid released from crystallizing melts derived from the lower crust or upper mantle. Intrusion of magmas below the exposed crustal level would also explain the high geothermal gradient calculated for the regional metamorphism.     

Fluid infiltration during contact metamorphism of interbedded marble and calc-silicate hornfels, Twin Lakes area, central Sierra Nevada, California

Abstract In the Twin Lakes area, central Sierra Nevada, California, most contact metamorphosed marbles contain calcite + dolomite + forsterite ± diopside ± phlogopite ± tremolite, and most calc-silicate hornfelses contain calcite + diopside + wollastonite + quartz ± anorthite ± K-feldspar ± grossular ± titanite. Mineral-fluid equilibria involving calcite + dolomite + tremolite + diopside + forsterite in two marble samples and wollastonite + anorthite + quartz + grossular in three hornfels samples record P± 3 kbar and T± 630° C. Various isobaric univariant assemblages record CO₂-H₂O fluid compositions of χCO₂= 0.61–0.74 in the marbles and χCO₂= 0.11 in the hornfelses. Assuming a siliceous dolomitic limestone protolith consisting of dolomite + quartz ° Calcite ± K-feldspar ± muscovite ± rutile, all plausible prograde reaction pathways were deduced for marble and hornfels on isobaric T-XCO₂ diagrams in the model system K₂O-CaO-MgO-Al₂O₃-SiO₂-H₂O-CO₂. Progress of the prograde reactions was estimated from measured modes and mass-balance calculations. Time-integrated fluxes of reactive fluid which infiltrated samples were computed for a temperature gradient of 150 °C/km along the fluid flow path, calculated fluid compositions, and estimated reaction progress using the mass-continuity equation. Marbles and hornfelses record values in the range 0.1–3.6 × 10⁴ cm³/cm² and 4.8–12.9 × 10⁴ cm³/cm², respectively. For an estimated duration of metamorphism of 10⁵ years, average in situ metamorphic rock permeabilities, calculated from Darcy's Law, are 0.1–8 × 10^?6 D in the marbles and 10–27 × 10^?6 D in the hornfelses. Reactive metamorphic fluids flowed up-temperature, and were preferentially channellized in hornfelses relative to the marbles. These results appear to give a general characterization of hydrothermal activity during contact metamorphism of small pendants and screens (dimensions ± 1 km or less) associated with emplacement of the Sierra Nevada batholith.     

海拉尔盆地乌尔逊凹陷古流体动力场数值模拟

首先对地下水动力场进行数学描述,以水动力为主线,建立适合该地区的数学模型,选择正确的地质模型参数,主要采用数值计算和模拟等手段,在恢复出乌尔逊凹陷主要目的层在不同地质历史时期古水动力的基础上又恢复出古流体势的演化历史。整体上,在乌尔逊凹陷,各主要地质时期的油势演化特征为西部地区大于东部地区,乌南大于乌北,且低势区的发育有较好的继承性。最后,依据油气聚集理论,进行油气运移聚集系统的划分、油气运移方向的分析,从而得出有关研究区油气运移聚集条件的认识,主要从流体动力的角度考虑,指出有利聚集区,认为在乌尔逊凹陷苏4～苏18井区、苏132～苏17～铜2井区及巴1～巴2～巴4井区为有利油气运移聚集区。     

湖北徐家山锑矿床流体包裹体特征及其意义

文章利用激光拉曼光谱和显微测温学方法,对湖北徐家山锑矿床成矿期的石英、重晶石和方解石中的流体包裹体进行了研究。研究表明,这些矿物中的流体包裹体主要有纯液体包裹体和液体包裹体(气相 液相)2类,其液体包裹体的气相成分为H2O±CO2±N2;石英、重晶石和方解石的均一温度分别为134~258℃、154~259℃和145~230℃,主要集中于150~200℃;流体的盐度w(NaCleq)和密度分别集中于3%~6%和0.90~0.96g/cm3。流体包裹体资料揭示出该矿床为典型的中低温热液锑矿床,其成矿流体为中低温、低盐度、中等密度热液。结合H、O、Sr、Pb同位素等研究结果,进一步推断该成矿热液主要是经深部循环演化的大气降水。     

新疆阿尔泰蒙库铁矿床的成矿流体及成矿作用

蒙库大型铁矿床赋存于上志留统—下泥盆统康布铁堡组变质火山-沉积岩系中,容矿岩石为石榴子石矽卡岩、变粒岩、浅粒岩和大理岩。矿体总体顺层分布,空间上与矽卡岩密切相关。研究表明,矽卡岩期石榴子石以发育玻璃质熔融包裹体、流体熔融包裹体和流体包裹体为特征,晚期矽卡岩阶段矿物中发育液相包裹体,变质期矿物中主要发育液相包裹体和含子矿物包裹体。矽卡岩期熔融包裹体的均一温度为1100℃,早期矽卡岩阶段流体包裹体均一温度变化于193～499℃,在450℃、350℃和230℃出现峰值。中期矽卡岩阶段均一温度变化于236～550℃,峰值为350℃。区域变质期均一温度介于132～513℃,在350℃、230℃和190℃出现峰值。流体包裹体的盐度w(NaCleq)介于1.23%～60.31%,流体密度变化于0.60～1.16g/cm3。石榴子石、石英和方解石的δ18OSMOW变化于0.2‰～8.4‰,δ18OH2O介于-5.1‰～5.33‰,δD为-127‰～-81‰,表明矽卡岩期成矿流体主要是岩浆水,混合少量大气降水;变质期流体主要为大气降水,为混合变质水。方解石δ13CPDB变化于-6.1‰～-2.3‰,表明流体中碳来自深部或地幔。成矿时代为早泥盆世早期(略晚于404～400Ma),成矿作用与矽卡岩的退化变质作用有关。     

豫西前河构造蚀变岩型金矿成矿过程中的流体-岩石反应

前河金矿区位于华北地台南缘,赋存在中元古界熊耳群安山岩和英安岩中,矿体受断裂破碎带控制。含矿热液在迁移过程中与围岩发生了广泛的流体-岩石反应而引起热液蚀变。本区石英中有4种类型的流体包裹体,均一温度范围为145～331℃,其中含CO2包裹体的完全均一温度主要分布在中-高温区。成矿流体的密度和压力变化范围分别是0.68～0.94g/cm3和(367.01～896.55)×105Pa。金大量沉淀成矿时的流体特征为:温度213～260℃、密度0.80～0.89g/cm3和压力(502.86～710.57)×105Pa。流体相为SO42->Na >Cl->K 型,CO2/H2O比值及N2、H2S、Ar、C2H6等挥发分的含量明显增高,f(CO2)、f(H2S)、f(CH4)和Eh值增大;f(O2)、f(H2O)和pH值减小。在青磐岩化安山岩的基础上发生的流体-岩石反应是造成本矿床金沉淀成矿的最主要原因。